Anisotropic photoconductivity in graphene 
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We investigate the photoconductivity of graphene within the relaxation time approximation. In 
presence of the inter-band transitions induced by the linearly polarized light the photoconductivity 
turns out to be highly anisotropic due to the pseudospin selection rule for Dirac-like carriers. The 
effect can be observed in clean undoped graphene samples and be utilized for light polarization 
detection. 
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INTRODUCTION 

Graphene membranes are optically transparent [l| 
as well as highly conductive^ even at room 
temperatures These two properties being incompat- 
ible with each other in conventional materials occur 
in carbon monolayers quite naturally and make them 
very promising for optoelectronical applications. 0, [B| 
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FIG. 1: Graphene Hall bar sample irradiated by linearly po- 
larized electromagnetic wave described by vector potential 
A. Applying a bias voltage leads to an electrical current j 
which depends on the photo-induced carrier concentration. 
The pseudospin orientation of the charge carriers described 
by Dirac Hamiltonian with the cone-shaped dispersion law 
shown by arrows is entangled with the particle momentum. 
The electrons in the valence band absorbing the photon en- 
ergy hv are excited to the conduction band producing the 
photoconductivity response. The electron-hole excitation rate 
is zero if the light is polarized along the pseudospins of the 
excited particles. In contrast, the excitation rate is maximal 
if the vector potential and pseudospin are perpendicular to 
each other. Since the pseudospin orientation is coupled with 
the particle's momentum the resulting photoconductivity er p h 
depends on the angle between A and j as shown in the inset. 
The absolute value of a ph is estimated in fig. [2] 



There is, however, another unusual property of carriers 
in graphene which makes this material even more 
interesting for optoelectronics. The carriers in graphene 
display an additional degree of freedom which is often 
dubbed as the pseudospin but, in fact, is connected to 
the sublattice index and has nothing to do with the real 
spin. @ We show, that the pseudospin manifests itself in 
the inter-band optical absorption making the transition 
probability sensitive to the pseudospin orientations in 
the initial and final states in a way similar to the real spin 
selective rules for the inter-band optical transitions in 
III-V semiconductors. Since the pseudospin is textured 
in the momentum space, as shown in fig. [TJ graphene's 
photoconductivity turns out to be anisotropic in the 
case of the linearly polarised light. The effect seems to 
be strong enough to find some applications in graphene 
optoelectronics. 



The model described below involves the optical exci- 
tation of the valence electrons to the conduction band 
of intrinsic (i.e. undoped) graphene. The idea is that 
the effective Hamiltonian describing the interaction be- 
tween the electromagnetic wave and carriers in graphene 
inherits the pseudospin-momentum entangled structure 
from the low energy kinetic term derived within the tight- 
binding approach. [6| Assuming normal incidence of a lin- 
ear polarized electromagnetic wave one deduces an elec- 
tron generation rate which strongly depends on the rel- 
ative orientation between the electron momentum and 
the linear polarization plane, see fig. [TJ As consequence, 
the photoconductivity is predicted to be anisotropic re- 
sulting in a high on/off ratio as a function of the lin- 
ear polarization angle. We note that the photoconduc- 
tivity in graphene has been also theoretically investi- 
gated in recent works, not analyzing its anisotropy. 
Moreover, the photoconductivity studied in this work 
should not be confused with the photocurrents predicted 
and measured [l3l - fl6j | in graphene. The pho- 
tocurrent can be generated without bias voltage applied, 
whereas the bias is necessary for the photoconductivity 
measurements. The photoconductivity and photocurrent 
anisotropy has been also found in the materials [l7l - fl9| 
other than graphene. 
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PRELIMINARIES 



The two-band effective Hamiltonian for 7r-system of 
graphene near half filling is Hq = hvp(a x k x + a y k y ), 
where vf ~ 10 6 ms _1 , k is the electron momentum, and 
o~ x>y are the Pauli matrices. The Pauli operator a rep- 
resents the pseudospin orientation which is depicted in 
fig. Q] for the eigenstates of H given by ^^{x.y) = 
^gik^x+ikyij ^ Ke je ) T , where tan6> = k y /k x , and k = ± 
denotes the band index, and the energy spectrum of Hq 
is E Kk = Khvpk. 

The interaction between the electromagnetic wave and 
charge carriers is described by the Hamiltonian H ln t — 
^-{o x A x + UyAy) and resembles the pseudospin struc- 
ture. Assuming the normal incidence and linear polar- 
ization of the electromagnetic wave A = Ao exp(— iut + 
ik z z) the golden-rule inter-band transition rate reads 



d 2 k'L 2 
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■wfVk', Kk){f K > k ' ~ Ink), (1) 



where f Kk is the distribution function, and 

x [5{E K , k , - E Kk -hw) + 5{E K , k , - E Kk + hw)} 
1 + kk' cos(0 + 9' - 26> pol ) 



(2) 



is the transition probability. Here u = 2itv is the ra- 
diation frequency, and tan# po i = A y /A x is the linear 
polarization angle. The length L plays a role of the sam- 
ple size or the laser spot diameter whichever is smaller. 
Eq. ([2]) describes the direct inter-band transitions and, 
thanks to the momentum and energy conservation, nat- 
urally includes ^-functions in the first two lines. Most 
important, however, is the third line which depends on 
the difference between the linear polarization angle 9 po \ 
and direction of carrier motion. This dependency disap- 
pears in the case of the circular polarization and is crucial 
for the effect considered below. 



PHOTOCONDUCTIVITY WITHIN THE 
RELAXATION TIME APPROXIMATION 



In the following we focus on the electron transport, 
i.e. K = +, and the carriers are excited from the valence 
to conduction band, as shown in fig. Q] To describe the 
recombination process we introduce the inelastic relax- 
ation time Tj which corresponds to the life time of the 
optically excited states. The steady state distribution 
function /i^ is then obtained by balancing the genera- 



tion rate (JTJ) and the relaxation rate /tj and reads 
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We naturally assume that the initial state is the equilib- 
rium one described by the Fermi-Dirac distribution func- 
tion There is no electrical current in the steady 
state described by the distribution function ([3]). 

The momentum relaxation is assumed to be due to the 
elastic scattering of carriers on impurities. The average 
momentum fiAk which the electrons gain due to the ex- 
ternal electric field E can be estimated as fiAk = eEr e , 
where r e is the elastic momentum relaxation time. For 
small electric field (linear response) the non-equilibrium 

(2) 

term fW can be obtained by expanding the steady-state 
function /jl^u with respect to small Ak in up to lin- 
ear order in E. Recalling hv = — <9Ak-E+(k-Ak) |Ak=o> the 
non-equilibrium distribution function for photo-excited 
electrons /V^ can be written as 
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Eq. Q is valid if and only if n 3> T e , i. e. optically 
excited states live much longer than the average time 
between two subsequent elastic scattering events. This is 
actually the case in graphene. 0, @ 

The current density due to the photo-excited electrons 
can be written as j p h = e J f^rv/T^ . This integral can 
be calculated in polar coordinates with the subsequent 
substitution e = E +k and reads 
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The photoconductivity for a given valley /spin channel is 
then given by 



. / 2 - cos(26> P oi) -sm(20 po i) , 
A P h ^ -sin(20 pol ) 2 + cos(20 po i) J W 



with the amplitude A p h being 



A ph = 



I6h 3 
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Rigorous analysis based on the Boltzmann equation writ- 
ten within the relaxation time approximation suggests 
the same expression for er p h but both t% and r e must be 
substituted by the total relaxation time r _1 = t~ 1 + 
t7 x + .... The effect of anisotropy predicted here does 
not depend on r anyway. Indeed, diagonalizing the ma- 
trix (j6]), the photoconductivity c ph = A p h parallel to the 
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light polarization plane turns out to be 3 times smaller 
than the perpendicular one a^ h = 3A p h, i.e. the pho- 
toconductivity is highly anisotropic, but the anisotropy 
itself is independent of r's. Thus, changing the linear po- 
larization angle from to 2tt one can observe two minima 
(and two maxima) in the current flow, as depicted in the 
inset of fig. [TJ These double extrema are a key signature 
of the effect predicted. 



DISCUSSION AND CONCLUSION 

Let us discuss the conditions necessary to observe the 
anisotropic photoconductivity given by eq. (j6|) and shown 
in fig. [TJ As it is clear from the analysis given in the previ- 
ous section the relative anisotropy does not depend on the 
relaxation times because the relaxation processes reduce 
the overall photoconductivity, not only its anisotropic 
part. The physical reason why the anisotropy does not 
vanish due to the momentum relaxation is the very fact 
that the anisotropic non-equilibrium distribution relaxes 
as fast as its isotropic contribution does. We believe 
therefore that the anisotropy can be detected easily as 
long as the photoconductivity response is large itself. 

To observe the photoconductivity the chemical poten- 
tial \i in graphene should be smaller than one-half of the 
excitation energy hjj/2 enabling direct excitations from 
the valence band. Assuming THz radiation, as used in 
the work by Karch et ai, [201 ] we arrive at the maximum 
\i less than lOmeV. Thus, the unintentional doping in 
graphene samples used before[i(| should be reduced by 
almost of two orders of magnitude. The temperature 
can also affect the effect even if the sample is perfectly 
neutral by reducing the photoconductivity by a factor 
of the order of hui/2T at zero chemical potential. Thus, 
room temperature T = 25meV seems to be somewhat 
to high for observing a sufficient signal at a radiation 
frequency of 1 THz. Moreover, the relaxation times r e 
and Tj assumed to be constant so far, will in fact also be 
temperature-dependent. However, one can facilitate the 
measurement by increasing the overall multiplier propor- 
tional to the radiation power, possibly by means of a high 
power pulsed NH 3 laser, plj 

In contrast to the photocurrents due to photon drag [II 



16, 20] the above effect is due to the pseudospin-selective 
inter-band transitions. The momentum transfer from 
photons to carriers is not important, and the effect should 
be observable even at normal incidence of light. The pre- 
dicted anisotropy is strongest for linearly polarized light 
source, whereas for circular polarization the transition 
probability ([2]) does not depend on the direction of car- 
rier motion, and the photoconductivity anisotropy does 
not occur. An elliptically polarized light source interpo- 
lates between these extreme cases. Moreover, the van- 
ishing anisotropy in the case of circular polarization can 
be used to separate the effect in question from the other 



photocurrent contributions. 1^- 16 , 2(| 

It is also interesting that the anisotropy predicted [12 1 
and observed [22[ recently in the photocurrent through 
graphene pn-j unctions seems to have the same origin as 
the one predicted here. There is, however, ir/2 off-set 
in the photocurrent vs. polarisation angle dependency 
as compared with the one shown in fig. [T] This is prob- 
ably because "the resulting photocurrent comes mainly 
from electrons moving nearly parallel to the barrier" 12| , 
and in order to maximize the concentration of such elec- 
trons the polarization plane must be set perpendicular to 
the pn-junction, i.e. along the photocurrent flow. One 
can reproduce this 7r/2 off-set also within our model by 
taking into account the dependence on k of the angle 9 
appearing in eq. (3) in the driving term of the Boltzmann 
equation. 

As already stated, the eigenvalues of the photoconduc- 
tivity tensor are predicted to differ by a factor of 3. In 
order to estimate the overall magnitude of the photocon- 
ductivity compared to other conduction mechanisms, let 
us compare the residual carrier concentration due to the 
unintentional doping with the one induced by the inter- 
band excitation. The former varies from 10 cm -2 for 
low mobility flakes on SiC>2 to 10 8 cm~ 2 for suspended 
samples after annealing. [23] The latter can be estimated 
as riph = Ti/(I/ 2 r p h) where r p h relates to the total photo- 
excitation rate as l/r p h = / d 2 ^2 ^[/+k]- On the other 
hand hoj / r p h can also be seen as the radiation energy ab- 
sorption rate which is nothing else than the absorbed ra- 
diation power W a . Note, that W a relates to the incident 
radiation power Wi as W a /Wi — na (where a — e 2 /(he) 
is the fine structure constant) for a single layer graphene 
membrane. [H [3] Thus, n p h at finite temperature T^O 
can be estimated as 



n ooq WiTi ft(0) 

nph = - 023 I 2 7^ (/ - 



(0)s 



(8) 



To be specific we assume that the photoconductivity is 
generated by a CH3OH laser (20]] with wavelength 118 jitm 
(i.e. huj = 10.5 meV) and Wi ~ 20 mW, and the sample 
itself is a suspended graphene membrane of the macro- 
scopic size slightly larger than the laser spot diameter 
of about 1 mm. Assuming r t ~ 1 ps 0, H[ we arrive at 
n ph - 2 • 10 7 cm~ 2 for L 2 ~ 1mm 2 and T = 0. This 
values are comparable to the residual carrier concentra- 
tion for suspended samples, 23] thus, the conductivity 
change in the irradiated graphene should be observable. 
Note, that n p h can be substantially increased by utilizing 
smaller samples and focusing the laser beam to a smaller 
spot. This requires a smaller radiation wave length (i. e. 
a higher laser frequency) to avoid diffraction effects. The 
results are summarized in fig. [2] 

The effect proposed above relies on the pseudospin tex- 
ture shown in fig. [TJ This texture remains stable as long 
as the low energy one-particle Hamiltonian Hq holds. At 
least from a theoretical point of view, the pseudospin 
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FIG. 2: The concentration of photoinduced carriers n p h as a 
function of temperature at different radiation frequency and 
sample size. (The latter is assumed to be roughly equal to 
the laser spot diameter and, therefore, has to be substan- 
tially larger than the radiation wavelength to avoid diffrac- 
tion effects.) The photoconductivity can be estimated as 
o"ph = e/in p h, where /j, is the mobility of carriers. The in- 
cident radiation power Wi and relaxation time Tj are 10 mW 
and 1 ps respectively. 



texture can be altered by electron-electron interactions 
which may be important in extremely clean samples. [25| 
This is the only fundamental obstacle for the photocon- 
ductivity anisotropy observation which we can see so far. 

To conclude, we predict strong anisotropy of the pho- 
toconductivity in graphene is presence of the linearly po- 
larized light. To observe the effect, we suggest to use 
undoped suspended graphene samples which allow the 
laser beam to excite the substantial number of photo- 
carriers from the valence band. The cleaner samples are 
expected to demonstrate the better results. They can be 
used as transparent detectors for the polarisation of the 
light passing through. 
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Echtermeyer for stimulating discussions. This work was 
supported by DFG via GRK 1570 and SFB 689. 



[1] R. R. Nair, P. Blake, A. N. Grigorenko, K. S. Novoselov, 
T. J. Booth, T. Stauber, N. M. R. Peres, and A. K. Geim, 



Science 320, 1308 (2008). 

[2] K. Bolotin, K. Sikes, Z. Jiang, M. Klima, G. Fudenberg, 
J. Hone, P. Kim, and H. Stormer, Solid State Communi- 
cations 146, 351 (2008), ISSN 0038-1098. 

[3] K. I. Bolotin, K. J. Sikes, J. Hone, H. L. Stormer, and 
P. Kim, Phys. Rev. Lett. 101, 096802 (2008). 

[4] P. Avouris, Nano Letters 10, 4285 (2010). 

[5] F. Bonaccorso, Z. Sun, T. Hasan, and A. C. Ferrari, Nat. 
Photon. 4, 611 (2010). 

[6] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 
(2009). 

[7] F. T. Vasko and V. Ryzhii, Phys. Rev. B 77, 195433 
(2008). 

[8] P. N. Romanets and F. T. Vasko, Phys. Rev. B 81, 
085421 (2010). 

[9] S. V. Syzranov, M. V. Fistul, and K. B. Efetov, Phys. 

Rev. B 78, 045407 (2008). 
[10] T. Oka and H. Aoki, Phys. Rev. B 79, 081406 (2009). 
[11] M. V. Entin, L. I. Magarill, and D. L. Shepelyansky, 

Phys. Rev. B 81, 165441 (2010). 
[12] S. Mai, S. V. Syzranov, and K. B. Efetov, Phys. Rev. B 

83, 033402 (2011). 
[13] F. Xia, T. Mueller, Y.-m. Lin, A. Valdes-Garcia, and 

P. Avouris, Nat Nano 4, 839 (2009). 
[14] X. Xu, N. M. Gabor, J. S. Alden, A. M. van der Zande, 

and P. L. McEuen, Nano Letters 10, 562 (2010). 
[15] J. Park, Y. H. Ahn, and C. Ruiz-Vargas, Nano Letters 

9, 1742 (2009). 

[16] J. Karch, P. Olbrich, M. Schmalzbauer, C. Zoth, C. Brin- 
steiner, M. Fehrenbacher, U. Wurstbauer, M. M. Glazov, 
S. A. Tarasenko, E. L. Ivchenko, et al., Phys. Rev. Lett. 
105, 227402 (2010). 

[17] S. K. Esayan, E. L. Ivchenko, V. V. Lemanov, and A. Y. 
Maksimov, JETP Lett. 40, 1290 (1984). 

[18] M. I. Karaman, V. P. Mushinskii, and G. M. Shmelev, 
Sov. Phys. Tech. Phys. 9, 730 (1983). 

[19] Y. S. Gal'perin and S. M. Kogan, Sov. Phys. JETP 29, 
196 (1969). 

[20] J. Karch, P. Olbrich, M. Schmalzbauer, C. Brinsteiner, 
U. Wurstbauer, M. Glazov, S. Tarasenko, E. Ivchenko, 
D. Weiss, J. Eroms, et al., Photon helicity driven electric 
currents in graphene, preprint arxiv: 1002. 1047. 

[21] S. Ganichev and W. Prettl, Intense terahertz excitation 
of semiconductors (Oxford University Press, 2006) . 

[22] T. J. Echtermeyer, priv. comm. 

[23] N. M. R. Peres, Rev. Mod. Phys. 82, 2673 (2010). 

[24] A. B. Kuzmenko, E. van Heumen, F. Carbone, and 
D. van der Marel, Phys. Rev. Lett. 100, 117401 (2008). 

[25] M. Trushin and J. Schliemann, Phys. Rev. Lett. 107, 
156801 (2011). 



